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ABSTRACT 

We present numerical simulations of an isothermal turbulent gas undergoing gravita- 
tional collapse, aimed at testing for "logatropic" behavior of the form Pt ~ log p, where 
Pt is the "turbulent pressure" and p is the density. To this end, we monitor the evolution 
of the turbulent velocity dispersion a as the density increases during the collapse. A 
logatropic behavior would require that a oc p~^^^, a result which, however, is not veri- 
fied in the simulations. Instead, the velocity dispersion increases with density, implying 
a polytropic behavior of Pf. This behavior is found both in purely hydrodynamic as 
well as hydromagnetic runs. For purely hydrodynamic and rapidly-collapsing magnetic 
cases, the velocity dispersion increases roughly as cr oc p^^^, implying Pt ~ p^, where Pt 
is the turbulent pressure. For slowly-collapsing magnetic cases the behavior is close to 
a oc p^/^, which implies Pt ~ p^^^. We thus suggest that the logatropic "equation of 
state" may represent only the statistically most probable state of an ensemble of clouds 
in equilibrium between self-gravity and kinetic support, but does not adequately rep- 
resent the behavior of the "turbulent pressure" within a cloud undergoing a dynamic 
compression due to gravitational collapse. Finally, we discuss the importance of the 
underlying physical model for the clouds (in equilibrium vs. dynamic) on the results 
obtained. 

Subject headings: ISM: clouds - magnetohydrodynamics - turbulence 
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1. Introduction 

Molecular clouds and clumps exhibit the well- 
known velocity dispersion- (or linewidth-)size relation 

a ^ R^/\ (1) 

where a is the linewidth-detcrmincd velocity disper- 
sion, and R the characteristic size. This correlation 
is observed both in ensembles of clouds (Larson 1981; 
Leung et al. 1982; Torrelles et al. 1983; Dame et al. 
1986; Myers & Goodman 1988; Falgarone, et al. 1992; 
Miesch & Bally 1994) or as a function of radius in 
quiescent cores using various tracers (Fuller & Myers 
1992; Caselh & Myers 1995; Goodman et al. 1997), 
although the latter studies have suggested that the 
scaling exponent in relation (|^) may actually differ 
between massive and low-mass cores. Furthermore, 
Goodman et al. have suggested that the exponent 
may decrease and approach zero as the innermost re- 
gions of the cores are considered, in which the turbu- 
lent velocity dispersion becomes subsonic. 

A second scaling relation, between mean density 
(p) and size, is also generally reported, reading 

{p)-R-\ (2) 

although its authenticity has been questioned on 
theoretical (Kegel 1989; Scalo 1990) and numerical 
(Vazquez-Semadeni, Ballesteros-Paredes & Rodriguez 
1997) grounds, and significantly discrepant scahng 
exponents have been found (e.g., Carr 1987; Loren 
1989), or none at all (e.g.. Plume et al. 1997). Rela- 
tions (H) and (H) constitute the now famous "Larson's 
relations" . 

In spite of the anomalies at small scales and high- 
mass regions, Larson's relations are generally ac- 
cepted as distinctive signatures of turbulence in molec- 
ular clouds and clumps (e.g., Larson 1981; Scalo 
1987), and together they imply 

(3) 

Relation (^) is actually a manifestation of virial equi- 
librium between the turbulent velocity dispersion (pos- 
sibly magnetohydrodynamic, or MHD) in the clouds 
(Larson 1981). A turbulent "pressure" Pt correspond- 
ing to the turbulent velocity dispersion can then be 
defined by (Lizano & Shu 1989; hereafter LS) 

= {dPt/dp), (4) 



in analogy with the relation between thermal pressure 
and the sound speed. Choosing 

Pt oc log p (5) 

recovers the virial relation (LS). Relation ^ is 
commonly refereed to as a "logatropic equation of 
state" 1^, or simply, a "logatrope" . 

It is important to emphasize that the concept of 
a turbulent "pressure" may not be a very realistic 
representation of the effects of turbulence, since it 
implicitly assumes a microscopic and isotropic pro- 
cess. Instead, turbulence is a phenomenon involv- 
ing a wide range of spatial scales, from the scale 
size of the system under consideration to the small- 
est dissipative scales. In particular, the existence of 
large-scale modes implies coherent motions which are 
more akin to ram pressure (locally anisotropic, with a 
well-defined direction) than to an isotropic, thermo- 
dynamic pressure. For these reasons, in the present 
paper we will focus primarily on the turbulent veloc- 
ity dispersion. References to the turbulent "pressure" 
will be made assuming that it can be defined accord- 
ing to eq. (^), for compatibility with published work, 
but the above caveat should always be kept in mind. 

The logatropic equation of state has been used in a 
number of studies of cloud support and stability, such 
as quasi-static contraction (LS), nonlinear wave prop- 
agation (Adams & Fatuzzo 1993; Adams, Fatuzzo & 
Watkins 1994; Gehman et al. 1996), and gravitational 
stability (McLaughlin & Pudritz 1996), among oth- 
ers. Nevertheless, the logatropic equation remains a 
completely empirical assumption, and there is no di- 
rect evidence that turbulent pressure indeed behaves 
in this manner in fully dynamic situations. In fact, 
there are some lines of reasoning which suggest that 
it may not: 

1. Larson's relations are observed in either ensembles 
of relaxed clouds (e.g., Torrelles et al. 1983; Dame et 
al. 1986; Myers & Goodman 1988; Falgarone, et al. 
1992; Miesch & Bally 1994), or as a function of ra- 
dius in quiescent cores (Fuller & Myers 1992; Caselli 
& Myers 1995; Goodman et al. 1997), but there is no 
evidence that they hold in fully dynamical processes, 
such as gravitational collapse. Interestingly, clouds 
which are strongly perturbed also seem to not follow 

^Strictly speaking, this is not an equation of state, since it does 
not invlove all three thermodynamic variables. However, we 
will allow ourselves the terminology for consistency with com- 
mon nomenclature. 
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the Larson scaling relations (e.g., Loren 1989; Plume 
et al. 1997). In general, there is little sampling of fully 
out-of-equilibrium, dynamical processes and, in par- 
ticular, dynamical collapse has not yet been directly 
detected because dynamical velocities occur only at 
very small scales. 

2. Although it might be argued that an ensemble of 
molecular clouds provides a complete sample of vari- 
ous dynamical stages, in actuality most observations 
refer to clouds and clumps close to equilibrium. Thus, 
the clouds included in surveys such as Larson's (1989) 
constitute an ensemble of equilibrium states for clouds 
of different masses rather than an ensemble of evolu- 
tionary steps for a single cloud (of constant mass). 
That is, instead of representing a number of different 
states for the same cloud, they represent the same 
state for different clouds. 

3. Vazquez-Semadeni et al. (1997) have suggested 
that there may exist large numbers of low-column 
density clouds which do not satisfy the density-size re- 
lation, and that possibly only the highest-column den- 
sity clouds follow such a scaling relation. Thus, the 
logatropic equation of state is not expected to apply 
to such low-column density clouds, which are prob- 
ably not in self-gravitating equilibrium, but rather 
pressure- or ram pressure-confined. 

As a first attempt to decide on this matter, in 
this paper we present two-dimensional numerical sim- 
ulations of a turbulent, self-gravitating, magnetized, 
isothermal gas, aimed at testing the variation of the 
velocity dispersion as a cloud is compressed by self- 
gravity. A related calculation has been performed by 
Bonazzola et al. (1987), who used low-resolution sim- 
ulations to estimate the correlation between the non- 
linear advection term (related to the turbulent pres- 
sure) and the density gradient in a compressible tur- 
bulent flow. 

We emphasize that the simulations discussed in 
this paper are not presented as models of cloud cores 
and their observed linewidths, but only as numerical 
"experiments" designed to test the applicability of the 
logatropic equation of state. Furthermore, through- 
out this paper we will refer exclusively to the non- 
thermal part of the velocity dispersion. In contrast 
with the observational situation, where the separation 
between the thermal and non-thermal components is 
an issue (e.g., Fuller & Myers 1992), in the simula- 
tions this is a trivial task, since there is no confusion 
between the fluid velocity and the thermal velocity 
dispersion, the latter being directly represented by 



the temperature field. 

The outline of the present paper is as follows. In 
we describe the numerical model; in ^ we present 
the results, for both purely hydrodynamic and fully 
MHD cases, and in §^ we summarize and discuss our 
results. 

2. Numerical model 

We numerically solve the full MHD equations in 
two or three dimensions (2D and 3D, respectively) in 
the presence of self-gravity for an isothermal, gravita- 
tionally-unstable gas, using the pseudo-spectral code 
described in Vazquez-Semadeni, Passot & Pouquet 
(1996), although here we restrict ourselves to a scale- 
free, isothermal case. The equations read 

1^ + V • ipu) ^ f^V'p, (6) 

i(V X B) X B - i^8V^u + 
P 

i^2(v2u+ivV-u), (7) 

— = V X (u X B) - i/gV^B -I- r/V^B, (8) 

V^ = p-l, (9) 

and 

P = c^p, (10) 

where, as usual, p is the density, u is the fluid velocity, 
P is the thermal pressure, B is the magnetic induc- 
tion, and 4> is the gravitational potential. The nondi- 
mensionalization is the same as in Passot, Vazquez- 
Semadeni & Pouquet (1995), to which we refer the 
reader for more details of the numerical method. The 
units are po = {p) (mean density in the integration do- 
main). Mo — Cs (isothermal speed of sound), Lq (size 
of the integration domain = 27r), to — Lajua (sound 
crossing time for the integration box), and Bo (mag- 
netic field strength such that at p = po, ^'A = Uo = Cg, 
where va is the Alfven speed). The resulting nondi- 
mensional parameters are J = Lq/Lj, the number of 
Jeans lengths in the integration box and the Mach 
number of the velocity unit = Uo/cg. Due to our 
choice of units, = 1. 
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Since the method is pseudospectral, it uses peri- 
odic boundary conditions and has no numerical dissi- 
pation. Due to the latter, dissipative operators need 
to be included explicitly. We use a combination of 
"hyperviscosity" and standard second-order viscosity 
in the momentum and magnetic equations which al- 
lows confinement of dissipative effects to the smallest 
scales in the simulation while filtering oscillations in 
the vicinity of strong shocks (Passot & Pouquct 1988; 
see also the discussion by Vazquez-Semadeni et al. 
1996). We also use a small amount of diffusion in the 
continuity equation (^) which also helps the code to 
handle strong shocks. Finally, the Poisson equation 
is used in a form suitable for handling infinite or 
periodic media (Alecian & Leorat 1988). 

The initial conditions for the simulations have a 
smooth Gaussian density profile peaked at the cen- 
ter of the integration box, with pmax = 3.35po and 
a FWHM of 0.53Lo. The initial velocity field is tur- 
bulent with Gaussian fluctuations of rms amplitude 
Wrms = 0.8cs and random phases, with exclusively ro- 
tational modes (i.e., no compressible motions). 

We have performed three non-magnetic simula- 
tions in two dimensions at resolutions of 128, 256 
and 512 grid points per dimension, respectively la- 
beled R128, R256 and R512. Except for the resolu- 
tion and the amount of dissipation (smaller at higher 
resolution), these runs are otherwise identical. A non- 
turbulent run at resolution of 512 grid point per di- 
mension labeled NT512 was also performed in order to 
test the numerical noise due to the grid discreteness in 
the calculation of the velocity dispersion (see below). 
Finally, a three-dimensional run, labeled run 3D96, 
with a resolution of 96 grid points per dimension was 
also performed in order to test for dimensional effects. 

Additionally, we performed four magnetic simula- 
tions. For these runs, labeled M256, MM256, M512 
and MM512, the initial magnetic field is along the a;- 
direction, with strength B^- = 0.2 for the "M" runs 
and Bx = 1.0 for the "MM" runs. For all runs except 
MM256 and MM512, the Jeans length is ij = 0.9Lo, 
while for MM256 and MM512, Lj = 0.67Lo. A sum- 
mary of the runs and their parameters is given in Ta- 
ble 1, including the diffusion coefficients. In all cases, 
r] = 0.002. 

The purpose of performing the three non-magnetic 
simulations at different resolutions is to test for con- 
vergence, in particular with regard to the effect of dis- 
sipation. Concerning the magnetic simulations, runs 
M256 and M512 are respectively similar to R256 and 



R512 except for the inclusion of the magnetic field. 
However, the magnetic field strength used in the "M" 
runs is quite small, so that the magnetic field does 
not prevent the gravitational collapse. Runs MM256 
and MM512 have a larger magnetic field (implying 
an Alfven speed equal to the sound speed) , somewhat 
closer to molecular cloud conditions. In order to guar- 
antee that the MM runs still undergo gravitational 
collapse, a smaller Jeans length was used. 

In the simulations we define the "collapsing cloud" 
as a circular region within the simulation, centered at 
the peak of the density distribution (calculated each 
time), containing 30% of the total mass. For the tur- 
bulent runs this may not be the best approximation, 
since the cloud shape is not really circular (see fig. 
1). However, using a true Lagrangian cloud bound- 
ary would be much more numerically involved, and we 
feel that our definition still provides reasonably accu- 
rate results, since ultimately gravity overpowers the 
turbulence and the shapes do not differ substantially 
from circular. For the 3D run, the cloud was defined 
as the region containing 10% of the total mass. 

The computation of the velocity dispersion requires 
some special care in order to remove the bulk infall 
velocity. This is a necessary step, since by definition 
the velocity dispersion is the root mean square ve- 
locity fluctuation, i.e., ((u — (u))^). However, in the 
present case of a collapsing cloud, the mean velocity 
is a function of radius Thus, we use the following 
procedure. We flrst compute the average infall speed 
Ur as a function of radius, and then compute the ve- 
locity dispersion as cr = ((u(r) — itr(r)f) )^/^, where f 
is the unit vector in the radial direction, and the aver- 
age is taken over the whole cloud, but using r = |r| at 
each position. This procedure shows why our Carte- 
sian grid introduces noise: the "circular" paths along 
which Ur is computed are not a perfect circumference, 
but rather the best possible approximation to one on 
a Cartesian grid, and the grid points on the path are 
not all at exactly the same distance from the center. 
Thus, even in the non-turbulent case there will be 
a systematic velocity dispersion at every radius, due 

^In general, the mean velocity will always be a function of po- 
sition in the case of a gas mass undergoing a global volume 
change. The simplest example is that of a gas mass in a cubic 
container being compressed by a piston on one side. In the 
direction of compression, the mean flow velocity will be a func- 
tion of position, being zero at the fixed wall, and equal to the 
velocity of the piston at the side of the piston. If the gas is ad- 
ditionally turbulent, the turbulent motions will be superposed 
on top of this mean-flow velocity. 
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to the presence of a radial velocity gradient and to 
the "thickness" of the circumference, the error being 
largest at the smallest radii. In order to estimate the 
magnitude of the numerical noise, we computed the 
numerical velocity dispersion in the nonturbulent case 
as well. 

3. Results 

As a typical example, figures 1 and 2 show the 
density and velocity dispersion fields at times t = 0.9 
and t — 2.4 in non-dimensional units for run R5I2, 
repectively. The former corresponds to the time of 
minimum velocity dispersion, after shocks have dissi- 
pated the initial velocity dispersion to a more slowly- 
decaying level, but before compression has begun to 
enlarge it again (see discussion about fig. 3 below). 
The latter time is the final state of the simulation, 
after which it stops because the code cannot handle 
the very large gradients developed at the center of the 
cloud any more. 

Figure 3 shows the log of the velocity dispersion vs. 
the log of the mean density for all runs as they evolve. 
For all the turbulent cases, it is seen that the initial 
transients suffer significant dissipation through shocks 
until a less dissipative regime is reached, at which the 
velocity dispersion is at a minimum. Subsequently, 
the velocity dispersion tends to increase with mean 
density, although with significant fluctuations. This 
behavior is in sharp contrast with relation (||). For 
the purely hydrodynamic runs R128, R256 and R5I2, 
a trend towards longer periods (i.e., larger density 
ranges) of nearly power-law behavior is noticeable. 
Although convergence on the duration of the power- 
law epoch of the evolution may have not been reached 
yet at 512^ resolution, the slope does appear to be 
converging to a value of 1/2. Even if convergence has 
not been attained yet, the observed trend is towards 
steeper slopes at higher resolutions, so in any case the 
discrepancy with relation (||) appears robust. 

An important possible problem is that this result 
might be an effect of the two-dimensionality of the 
simulations. Run 3D96 was performed as an attempt 
to resolve this question, although the resolution is 
necessarily lower. The evolution of the velocity dis- 
persion and the mean density for this run is also 
shown in fig. 3. Although at much slower rates than 
in the 2D runs due to the higher dissipation inher- 
ent to the lower resolution, the trend in run 3D96 
is still towards increasing a with (p) after the initial 



transients have passed. Thus, even though run 3D96 
does not permit confirmation of the rates approached 
by the high- resolution 2D runs, the increasing trend 
of a with (p) is maintained, suggesting that this be- 
havior is real, rather than just an effect of the two- 
dimensionality. 

In this regard, note that in general the simulations 
overestimate the viscous dissipation rate, since, for 
numerical reasons, the viscous coefficients have to be 
chosen so that the dissipative scales fit within the res- 
olution of the simulation. Instead, in the actual in- 
terstellar gas, the dissipation scales may be many or- 
ders of magnitude smaller than the scales of interest. 
Thus, the dissipation rate of a found in the simula- 
tions is at worst a lower bound to the actual rate, 
and the net increase of a found in the simulations is 
expected to be a real effect. In particular, run 3D96 
is the most dissipative of the runs performed, but a 
net increase is found also in this case. 

The non-turbulent run NT512 also exhibits a ve- 
locity dispersion which increases with mean density. 
This should be interpreted as a numerically-generated 
velocity dispersion which increases at larger infall 
speeds because the radial velocity gradient is also 
larger. Nevertheless, this numerical velocity disper- 
sion is seen to be generally about two orders of mag- 
nitude smaller than that for the turbulent runs (note 
that the curve for the non-turbulent run has been dis- 
placed upwards by an amount of 1.5 in log cr so that 
it fits within the plot). Thus, we rule out numeri- 
cal noise as the cause for the trends observed for the 
turbulent runs in fig. 3. 

The magnetic runs also exhibit a trend of increas- 
ing velocity dispersion with increasing mean den- 
sity, although with quite stronger fluctuations. Runs 
MM256 and MM512 exhibit a range of densities for 
which again roug hly a - j^uns M256 and 

M512 exhibit a somewhat slower dependence, close 
to cr ^ p^/*. In any case, the general trend is the 
same as in the non-magnetic cases, contrary to the 
logatropic behavior, relation (|^). 

4. Conclusions 

4.1. Summary and discussion 

We have argued that Larson's (1981) relations and 
the resulting logatropic "equation of state" (relation 
[^) and virial condition (relation [^) may describe an 
ensemble of clouds in (near) equilibrium between self- 
gravity and the turbulent velocity dispersion, but not 
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out-of-equilibrium, dynamical processes occurring on 
a single cloud. We have tested this assertion by means 
of numerical simulations of collapsing clouds with an 
initial turbulent velocity field, in both magnetic and 
non-magnetic regimes. 

The simulations exhibit in all cases a turbulent 
velocity dispersion which increases with mean den- 
sity as the collapse proceeds, in contradiction with 
the expected behavior for a logatrope, relation (|^). 
Non- magnetic and strongly self-gravitating runs seem 
to approach a power-law behavior of the form a ~ 
(/o)^/^, while weakly self-gravitating magnetic runs in 
general tend to have shallower dependences, although 
always with positive exponents. In particular, the 
fact that magnetic runs exhibit the same qualitative 
behavior suggests that weak magnetic fields cannot 
induce a logatropic (or a cr behavior either. 

Interestingly, run M512 exhibits a behavior very close 
to fj ~ p^/^. Assuming that this run has converged to 
the true slope, it is noteworthy that the implied tur- 
bulent pressure satisfies Pt p'^^^ . This result is con- 
sistent with that of McKee & Zweibel (1995) for the 
polytropic index of Alfven waves under slow compres- 
sion. However, runs MM256 and MM512 appear to 
be closer to the a p^/^ (Pt oc p^) behavior observed 
in the non-magnetic runs. This distinction is likely to 
be due to the larger Jeans length used in the M runs, 
implying a slower collapse (final time tfin = 2.4) than 
for the MM runs (tfin = 1.0), so that the M runs are 
closer to the slow compression assumption of McKee 
and Zweibel. 

We emphasize that although convergence may not 
have been fully achieved yet at the highest resolu- 
tion we used (512 x 512 grid points), the trend is 
towards faster increase of the velocity dispersion with 
density at higher resolution, away from the behavior 
predicted by the logatropic equation. Thus, the re- 
sult that the velocity dispersion increases with mean 
density appears quite robust. Moreover, the trend of 
increasing a with {p) is preserved in run 3D96 (albeit 
at a slower rate due to the lower resolution of this 
run), thus ruling out the possibility that our results 
are a purely 2D effect. 

The main consequence of our results is that the 
logatropic "equation of state" appears to be inade- 
quate for the description of dynamical processes oc- 
curring in a cloud. This implies that the use of the lo- 
gatropic equation in studies of gravitational collapse 
and dynamical stability is questionable. Its use in 
studies of quasi-static processes (e.g., LS) may still 



be justified, although in general the question remains 
open as to whether the logatropic equation, which can 
be thought of as representing the final states of the 
virialization process, also represents the behavior of 
the turbulent pressure during the relaxation processes 
which lead to virialization. For this reason, it would 
also be interesting to test its applicability in problems 
of nonlinear wave propagation (e.g., Adams & Fatuzzo 
1993; Adams, Fatuzzo & Watkins 1994; Gehman et 
al. 1996). 

Finally, we remark that the ensemble consisting of 
the evolutionary states of our simulated gravitationally- 
collapsing clouds with a fixed mass is completely dif- 
ferent from the ensemble constructed from the ob- 
servations of many clouds of different masses in near 
equilibrium. The present work shows that for the for- 
mer ensemble, the logatropic equation of state is not 
applicable. 

4.2. Comparison vifith previous work 

In our simulations we obtain a polytropic form 
{Pt oc p'^'") for the effective "equation of state" of 
the turbulence, with polytropic exponents 7cff = 2 for 
the non-magnetic and strongly self-gravitating cases, 
and 7off = 3/2 for the weakly self-gravitating cases. 
This result appears to be in contradiction with that 
of McLaughHn & Pudritz (1996, hereafter MP), who 
conclude that the total pressure (thermal plus turbu- 
lent) is not expected to behave as a polytrope. MP 
reach this conclusion on the basis of a stability anal- 
ysis, noting that truncated polytropic solutions with 
< 7eff < 1 (consistent with the observed lower tem- 
peratures of denser structures) have never unstable, 
or even critically stable, equilibrium solutions. That 
is, absolutely stable configurations are discarded by 
MP so that a cloud is eventually able to collapse and 
form a star. 

It can then be seen that the difference between our 
results and those of MP arises from the considera- 
tion of different physical models for the clouds. While 
MP's clouds arc in hydrostatic equilibrium, our clouds 
are always out of equilibrium and are already unsta- 
ble from the start. These may originate from clumps 
rendered unstable by external turbulent compressions 
(Vazquez-Semadeni et al. 1996) if the effective equa- 
tion of state (i.e., the heating and cooling) permits it. 
In such cases, the clumps never need to pass through 
a static equilibrium state. Another possibility is the 
well-known onset of gravitational instability due to 
the loss of magnetic support caused by ambipolar dif- 
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fusion (e.g, Nakano 1979; LS). 

Finally, we note that the polytropic exponents im- 
plied by our simulations are larger than the critical 
value 7c =4/3 below which gravitational collapse can 
proceed to a singularity (e.g., Chandrasekhar 1961). 
Thus, if turbulent pressure continued with this be- 
havior unrestrictedly, it would eventually halt the col- 
lapse. However, we do not expect this occur since, at 
late stages of the collapse, dissipation becomes im- 
portant again due to the large velocity gradients that 
develop. In fact, in fig. 3 an end to the steady in- 
crease of a is seen at large values of (p) for several of 
the runs. Thus, we speculate that turbulent pressure 
cannot by itself halt the collapse. 
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Run 




FWHM/L^ ' 


Lj/L), 










(h) 




R128 


3.35 


0.53 


0.9 


0.8 





0.03 


2.00 X 10-3 


8.12 X 10" 


■12 


R256 


3.35 


0.53 


0.9 


0.8 





0.0075 


5.00 X 10-'' 


3.13 X 10" 


-14 


R512 


3.35 


0.53 


0.9 


0.8 





0.008 


1.25 X 10-* 


2.00 X 10" 


■16 


NT512 


3.35 


0.53 


0.9 








0.008 


1.25 X 10-^ 


2.00 X 10" 


16 


M256 


3.35 


0.53 


0.9 
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0.2 


0.0075 


5.00 X 10-4 


3.13 X 10" 
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3.35 


0.53 
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1.25 X 10-* 


2.00 X 10" 
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3.35 


0.53 


0.67 
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0.0075 


5.00 X 10-4 


3.13 X 10" 


■14 


A1M512 


3.35 


0.53 


0.67 


0.8 


1.0 


0.008 


1.25 X 10-4 


2.00 X 10" 


16 


3D96 


4.93 


0.54 


0.9 


0.8 





0.013 


1.30 X 10-3 


4.00 X 10" 
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(a) Central density of initial density peak in units of 
mean density 

(b) FWHM of initial density peak in units of box size 

(c) Jeans Icngtli in units of box size 

(d) Initial runs turbulent speed in units of sound 
speed 

(e) Initial strength of uniform magnetic field 

(f) Diffusion coefficient for the continuity equation 

(g) Standard viscosity coefficient 

(h) Hyperviscosity coefficient 



8 




Fig. 1. — Views of the density (top) and velocity dis- 
persion {bottom) fields for run R512 at time t = 0.9 in 
non-dimensional code units. This is the time of min- 
imum velocity dispersion, having been already dissi- 
pated by shocks and not yet enhanced by the gravita- 
tional compression. The velocity dispersion is shown 
within the circle contining 0.3 of the total mass in the 
simulation. The gray scale for the density is logarith- 
mic. 




Fig. 2. — Same as fig. 1 at time t = 2.4. This is the 
final stage of the collapse. Note that the gray scale 
in this figure differs (has a larger maximum value) 
from that in fig. 1 in order to maximize clarity. The 
maximum density is Pmax = 330 (p). 
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Fig. 3. — Log-log plot of the velocity dispersion vs. 
average cloud density for all runs. The plot for run 
NT512 has been displaced upwards by an amount of 
1.5 in the y-axis units. In all cases, the velocity dis- 
persion is seen to increase with average density. The 
circles in the curve for run R512 indicate the individ- 
ual times during the collapse at which the curve was 
sampled, taken at intervals At = 0.1 code time units. 
The solid straight line shows a power-law with expo- 
nent 1/2, and the dotted straight line shows a slope 
of 1/4. 
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